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ABSTRACT 

Iron emission lines at 6.4 — 6.97 keV, identified with Ka radiative transitions, 
are among the strongest discrete features in the X-ray band. These are one of 
the most powerful probes to infer the properties of the plasma in the innermost 
part of the accretion disk around a compact object. In this paper we present 
a recent Suzaku observation, 100-ks effective exposure, of the atoll source and 
X-ray burster 4U 1705-44, where we clearly detect signatures of a reflection 
component which is distorted by the high -velocity motion in the accretion 
disk. The reflection component consists of a broad iron line at about 6.4 keV 
and a Compton bump at high X-ray energies, around 20 keV. All these features 
are consistently fitted with a reflection model, and we find that in the hard 
state the smearing parameters are remarkably similar to those found in a 
previous XMM-Newton observation performed in the soft state. In particular, 
we find that the inner disk radius is Ri n = 17 ± 5 R g (where R g is the 
Gravitational radius, GM/c 2 ), the emissivity dependence from the disk radius 
is r _2 ' 5±0 ' 5 , the inclination angle with respect to the line of sight is i = 43°±5°, 
and the outer radius of the emitting region in the disk is R ou t > 200 R g . We 
note that the accretion disk does not appear to be truncated at large radii, 
although the source is in a hard state at ~ 3% of the Eddington luminosity for 
a neutron star. We also find evidence of a broad emission line at low energies, 
at 3.03 ± 0.03 keV, compatible with emission from mildly ionized Argon (Ar 


2 T. Di Salvo et al. 


XVI-XVII). Argon transitions are not included in the self-consistent reflection 
models that we used and we therefore added an extra component to our model 
to fit this feature. The low energy line appears compatible with being smeared 
by the same inner disk parameters found for the reflection component. 

Key words: line: formation line: identification — stars: neutron — stars: 
individual: 4U 1705-44 — stars: magnetic fields — X-ray: general — X-ray: 
binaries 


1 INTRODUCTION 


Neutron Star Low Mass X-ray Binaries (hereafter NS LMXBs) are binary systems in which 
a weakly magnetic NS accretes matter from a low mass (< 1 M & ) companion star via Roche- 
Lobe overflow. In these systems the accretion disk can approach the compa ct object, as 
testi fied by the very fast time variability observed up to kHz frequencies (see 


2006 


van der Klis 


, as a review). Broad emission lines (FWHM up to ~ 1 keV) at energies in the 


range 


Piraino et al. 

2007; 

Bhattacharvva & Stro 

imaver 

2007; 

Cackett et al. 

2008[; 

Pandel et al 


2008; 

di Salvo et al 

2009( 

D’Ai et al. 

2009; 

Iaria et 

a 

2009; 

Paoitl 

:o et a 

1. 2009 

Caeke 

t et al. 

2009; 

Shanos 

mikov 

et al. 

2009; 

Paoitto et al. 

201( 


daron et al. 

2011; 

Piraino et al. 

2012 



Papitto et al. 


2013 


Miller et al. 


20131. These lines are identified with the Ka radiative tran¬ 


sitions of iron at different ionization states. These features are powerful tools to investigate 
the structure of the accretion flow close to the central source; in particular, important in¬ 
formation can be obtained from the detailed spectroscopy of the line profile, since it is 
determined by the ionization state, geometry and velocity field of the reprocessing plasma. 
In fact, the broad iron line observed in NS LMXBs is thought to originate from reflection 
of the primary X-ray continuum off the inner accretion disk and the width of the line is 
expected to be a signature of the Keplerian motion of matter in the inner accretion disk at 
(mildly) relativistic velocities. In this model, the combination of Doppler effects from the 
high orbital velocities and Special and General relativistic effects arising from the strong 
gravity in the vicinity of the NS smears and shifts the reflected features. As a consequence, 
the line will have a characteristically broad and asymmetric profile, the detailed shape of 
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which depends on the inclination and on how deep th e acc retion disk extends into the NS 
gravitational potential (e.g. 


Fabian et al. 

1989; 

Matt et al. 

1992) 


If the origin of this line is from disk reprocessing, one would also expect the presence 
of other discrete features (such as emission lines and absorption edges from the most abun¬ 
dant elements) and a bump between 20 and 40 keV due to direct Compton scattering of 


observed in the spec 

tra of some IN 

S LM 

XBs (see e.g. 

Bar re 

et al. 

2000; 

Piraino et al. 

1999; 

Yoshida et al. 

1993; 

Fiocchi et al. 

2007; 

Egron et al. 

2013; 

Miller et al. 

2013 

), usually with 


reflection amplitudes (defined in terms of the solid angle Q/2tt subtended by the reflector 
as seen from the corona) lower than 0.3. In some cases an anti-correlation has been claimed 


between the photon ind ex of the primary spectrum and 


reprocessed component (Zdzi ar ski et al. 


1999 


Barret et al. 


ie_reflection amplitude of the 


2000 


Piraino et al. 


19 991 ). the 


same observed in Seyfert galaxies and galactic Black Hole (hereafter BH) candidates. This 
is probably caused by variations in the position of the inner rim of the disk. 

The disk origin of the iron line in NS LMXBs is, however, debated in literature because 
of the bright ness o f these sources, which may cause photon pile-up and systematics in CCD 
spectra (INg et al.ll2010 l). making the detection of any asymmetry in the line profile somewhat 
controversial. However, a large number of simulations using a statistical model of photon 
pile-up to assess its impacts on relativistic disk line and continuum spectra suggest that 
severe photon pile-up acts to falsely narrow emission lines, leading to falsely large disk radii 
(Miller et ah 2010J). These simulations also indicate that relativistic disk spectroscopy is 
generally robust against pile-up when this effect is modest. Moreover, several authors (e.g. 


Cackett et al. 


2012 


Kgron et al 


20131 ) have shown that the CCD-based spectra from Suzaku 


and XMM-Newton are compatible with gas-based spectra from EXOSAT, BeppoSAX, and 
RXTE, demonstrating that the broad profiles seen are intrinsic to the line and not due to 
instrumental issues. They also report that a few BeppoSAX observations show evidence for 
asymmetric lines, with a relativistic di skline model providing a significantly better fit than 


a Gaussian line (see also 


Piraino et al. 


1999). 


Nevertheless alternatives have been proposed to explain the profiles of these features. In 
particular, INg et al.l (20 101) propose that Compton broadening may be sufficient to explain 
the large width of the line. However, when self-consistently included in the fit, Compton 
broadening alone appears to be insufficient to expla in the observed line pro fi le (se e e.g. 


Reis et al. 


2009 


Eg ron e t a l. 20131: 


Sanna et al. 


20131) . Also, 


Cackett fc Miller 


(2013) 


have 
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explored the observational signatures expected from broadening in a wind. In this case the 
iron line width should increase with increasing the column density of the absorber (due to 
an increase in the number of scatterings). They show that there is no significant correlation 
between line width and column density, favoring an inner disk origin for the line broadening 
rather than scattering in a wind. 

4U 1705-44 is a well-studied atoll source (see 


Has inger fe van der Klisl 1989 ; 


Olive et al. 


2003]), which also shows type-I X-ray bursts. Similarly to X-ray binaries containing BHs, 


this source regularly shows state transitions: from a high/soft state, where the X-ray spec¬ 
trum is dominated by soft spectral components with typical temperatures less than a 
few keV, to low/har d states where the X-ray spectrum is do minated by a hard thermal 
Comptonisation (e.g. iBarret fe Olivel l2002t iPiraino et al.l l2007f) . The presence of broad dis¬ 


crete features in this source has been often reported in literature. A broad iron line was 


observed, in the soft and/or hard_s 


Chan dra/HETG ( Di Salvo et al. 2005), BeppoSAX (Piraino et al 


20 0 9 ). and XMM/pn ( Di Salvo et al 


ate, with mo derate/high spectr al resoluti on by the 


2007 ), Suzak u ( Reis et al. 


2009; 

D'Ai et al. 

2010; 

Eeron et ah 

2013) 


observation, taken in August 2008 during a soft state (45 ks effective exposure time), showed 
one of the highest signal-to-noise ratio iron line profile ever observed in a NS LMXB. The 
line profile is cle arly broad and cou ld be fitted equally well with a relativistic line profile, 


such as diskline (IFabian et al.lll989l) or relline (IDause 


r et ah 


reflection models, such as reflionx (IRoss fc Fabianll2005l ). refbb (IBallantvnell2004Ji . and xillver 


2010). or with self-consistent 


(iGarcfa & Kallman 2010). All these models gave parameters of the inner disk with unprece¬ 


dent ed precision and all compatible with each other within t 


see 


Di Salvo et al. 


2009 


D’Ai et al. 


20101 


Fgron et ah 


re small statistical uncertainties 


201 3). The line is identified with the 


K a transition of highly ionised iron, Fe XXV; the inner disk radius is Ri n = 14 ±2 R 9 (where 
R 9 is the Gravitational radius, GM/c 2 ), the emissivity index of the disk is —2.27±0.08 (com¬ 
patible with a disk illuminated by a central source), the inclination angle with respect to 
the line of sight is i = 39 ± 1 degrees. This, together with the presence of other low-energy 
features from S XVI, Ar XVIII, Ca XIX and a smeared iron edge at 8.4 keV, which all are 
compatible with being smeared with the same inner disk parameters, makes 4U 1705-44 
the best source for a detailed spectroscopic study, in order to address the disk origin of the 
observed iron line and of the whole reflection component. 

In this paper we present a high statistics, 100-ks effective exposure, Suzaku observation 
of 4U 1705-44 during a hard state: these data allow us a detailed study of the reflection 
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features and the fit, with a self-consistent reflection model, of both the iron line profile 
and the associated Compton reflection bump at energies above 10 keV. In this spectrum, 
which includes hard-band data (up to 200 keV), the overall fractional amount of reflection is 
well determined by fitting the Compton bump. We can therefore test whether the observed 
iron line is consistent with this fractional amount of reflection. In this way we confirm 
independently (fitting a different spectral state and using different instruments) the inner 
disk parameters already obtained with XMM-Newton in the soft state. 


2 OBSERVATIONS 


Suzaku ( Mitsuda et all 2007) observed 4U 1705-44 on 2012 March 27 as the result of a 


Target of Opportunity (ToO) program during a hard state for a total observing time of 250 
ks, corresponding to an effective exposure time of about 100 ks because of observational gaps 


caused by Earth occultations along the low equa' 
X-ray Imaging Spectrometers (XIS, 0.2-12 keV; 


Detector (HXD, 10-600 keV; 


Takahashi et ah 


orial orbit of the Suzaku satellite. Both the 


Kovama et ahl (120071 )) and the Hard X-ray 


(2007)) instruments were used during these 


observations. There are four XIS detectors, numbered as 0 to 3. XIS0, XIS2 and XIS3 all use 
front-illuminated CCDs and have very similar responses, while XIS1 uses a back-illuminated 
CCD. At the time of this observation the available CCDs were three due to the loss of the 
XIS2. The HXD instrument consists of two types of collimated (non-imaging) detectors, the 
PIN diodes (10-70 keV) and the GSO scintillators (30-600 keV). 

We reprocessed the data using the aepipeline tool provided by Suzaku FTOOLS version 
20 updated with the latest calibration hies (2013 November). As second step, in order to 
obtain a more accurate estimate of the Suzaku attitude, we have calculated a new attitude 
using the free tool aeattcor.sl created by J. E. Davis. Then we have applied the new attitude 
to XIS event hies using the FTOOLS xiscoord. During the observation, XIS0, XIS1, and 
XIS3 were operated using the 1/4 window option. The effective exposure time of each XIS 
CCD is 96.67 ks. In order to estimate the pile-up in the XIS spectra we have used the 
public available tool pile-estimate.si created by M. A. Novak. Using a circular region with 
radius equal to 105", we have found that in each XIS CCD, the pile-up fraction is ~ 3% at 
most. The pile-up fraction is sufficiently small that we can neglect its effects on our spectral 
fitting results. In fact, we have checked that spectral results do not change signihcantly if 
we exclude a central circle in the extraction region in order to further reduce the pile-up 
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fraction. Therefore, we have extracted the source and background spectra from a circular 
extraction region of radius 105" each, the background circle being centered close to the edge 
of the CCD, where no significant contaminating photons from the source were present. The 
response hies of each XIS spectrum have been generated using the xisrmfgen and xisarfgen 
tools. Since the response of XISO and XIS3 are very similar, we have combined their spectra 
and responses using the tool addascaspec. 

The PIN spectrum has been extracted using the tool hxdpinxbpi. Both the non X-ray 
and cosmic X-ray backgrounds are taken into account. The non X-ray background (NXB) is 



by the HXD team are used. We selected the HXD/PIN events in the energy range 12-30 keV 
and produced the HXD/PIN background-subtracted light curve using the SUZAKU tool 
hxd-pinxblc and adopting the background event hies distributed by the HXD team. Since 
this light curve follows the XISO light curve, we conclude that no contaminating hare was 
present in the data. We have also extracted the GSO spectrum running the tool hxdgsoxbpi. 
For the background we have used the ’tuned’ non X-ray background, whereas for the response 
hie we have used the latest version provided on 2011-06-01. 

We have extracted the XISO light curves in the energy range 0.9 — 2.8, 2.8 — 10 and 
0.9 — 10 keV, respectively (see Fig. 1). Nine type-I X-ray bursts are observed in the total 
250-ks light curve. The source (persistent) count rate gradually increases by about 40% 
during the observation. Since no changes are observed in the hardness ratio (given by the 
ratio of the source counts in the 2.8 — 10 keV range to the source counts in the 0.9 — 2.8 
keV, see Fig. 1) we conclude that the X-ray spectral shape of the source does not change 
significantly during the observation. 


3 SPECTRAL ANALYSIS AND RESULTS 

In order to extract spectra for the persistent emission, we have excluded the type-I bursts 
that occurred during the observation. In particular we have excluded approximately 100 s 
of data starting from the onset of each burst. We adopt 0.7 — 11 keV energy range for the 
XIS0+XIS3 (hereafter XIS03) and XIS1 spectra, 15 — 50 keV energy range for the HXD/PIN 
spectrum and 50 — 200 keV energy range for the HXD/GSO spectrum. We excluded the 
energy interval between 1.7 and 2.0 keV from the XIS03 and XIS1 spectra because of the 
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Figure 1. Left: Suzaku XISO light curve in the energy range 0.9 - 10 keV showing two of the nine type-I bursts which 
occurred during the 250-ks observation. Right: Suzaku XISO light curves in the energy range 0.9 — 2.8 keV (top panel), 
2.8 — 10 keV (middle panel), and the corresponding hardness ratio (bottom panel). 


presence of systematic features associated with neutral silicon and neutral gold which give a 
mismatch between the two spectra. The XIS spectra were grouped by a factor 4 in order not 
to oversample too much the instrumental energy resolution. The HXD/PIN and HXD/GSO 
spectra were grouped in order to have at least 25 photons per energy channel. We fitted the 
spectra using XSPEC version 12.7.0. 

We started to fit the continuum in the 0.7-200 keV energy range with the typical model 


used for NS LMXBs o 
this source too (see e.g. 


the atoll class, w' 


rich revealed 


Di Salvo et al. 


2009 


Piraino et al. 


o be the best fit continuum for 


2007|; 


Barret fc Olive 


2002). This 


model consists of a soft blackboc 


modelled by nthComp (iZvcki et ah 


v and a thermal Comptonised component, in this case 


19991 ), modified at low energy by photoelectric absorption 


caused by neutral matter and modeled by phabs in XSPEC. This continuum model gave, 
however, an unacceptable fit, corresponding to a x 2, /dof = 2425.13/1511, because of the 
presence of evident localised residuals at 2.5-3.5 keV, 6-9 keV and 15-30 keV. The most 
prominent is a clear iron line profile at energies from 5 to 7 keV and an absorption feature 
at 7-8 keV (see Fig. 2). 

In order to fit these re_siduals, we first added to our continuum model the pexriv com¬ 
ponent (Magdziarz & Zdziarskj 1995 ) which takes into account the iron edge and Compton 
bump present in the residuals. Note that pexriv does not self-consistently include any emis¬ 
sion line. The photon index and the high energy cutoff of the pexriv component were fixed 
to those of the nthComp component. Here and in the following we neglect any reflection 
of the soft (blackbody) component, which in any case contributes to a small fraction of 
the total flux and most of its flux is at soft energies (below 4 — 5 keV). The addition of 
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Energy (keV) 


Figure 2. Suzaku data in the energy range 0.7 - 200 keV (top) and residuals in units of a with respect to the simpler 
phenomenological model (bottom) of 4U 1705-44. The model consists of a blackbody (dotted lines) and the Comptonization 
component nthComp , both multiplied by photoelectric absorption. 


this component gave a significant improvement of the fit reducing the x 2 /dof to the value 
2009.67/1509 (Ay 2 = 415 for the addition of two parameters). Therefore, the presence of the 
Compton hump and the iron edge are detected with high statistical significance (an F-test 
would give a probability of chance improvement of the fit negligibly small, ~ 2.7 x 10 -62 ). 
We then added a Gaussian line at 6.4 keV obtaining again an improvement of the fit. As 
a first step we decided to fix the energy of the line at 6.4 keV because, otherwise, it tends 
to get broad (Gaussian sigma about 1 keV) and its energy tends to decrease to 6.2 keV. 
This fit gives a x 2 /dof = 1956.50/1507 (Ay 2 = 53 for the addition of two parameters). This 
line can be identified with fluorescence from neutral iron. The addition of another Gaussian 
at about 3 keV again improves the fit, giving y 2 /do/ = 1829.89/1504 (Ay 2 = 127 for the 
addition of three parameters). This line can be identified with emission from mildly ionized 
Argon, Ar XVI-XVII. In order to check whether the iron line energy was stable or not, 
we let the iron line energy free to vary obtaining y 2 jdof = 1818/1503. We also tried to 
substitute the Gaussian line at ~ 6 keV with a diskline. In this case we had to fix all the 
smearing parameters but the inner radius of the disk. This did result in a slight improvement 
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Table 1. The best fit parameters of the spectral fitting of the Suzaku spectrum of 4U 1705-44 in the 0.7 — 200 keV energy 
range with phenomenological models. The blackbody luminosity is given in units of L 35 /D 1 Q, where L 35 is the bolometric 
luminosity in units of 10 35 ergs/s and Dio the distance to the source in units of 10 kpc. The blackbody radius is calculated in 
the hypothesis of spherical emission and for a distance of 7.4 kpc. Fluxes in the nthComp and pexriv components are calculated 
in the 1 — 16 keV range, while total flux is calculated in the 1 — 10 keV band. Uncertainties are given at 90% confidence level. 


Component 

Parameter 

Basic Model 

Pexriv 

Pexiriv + 2Gauss 

Smeared (Pex + 2Gaus) 

phabs 

N h (xio 22 cm -2 ) 

1.899 ±0.025 

2.13 ±0.05 

2.02 ±0.05 

2.04 ± 0.06 

bbody 

kT BB (keV) 

0.503 ± 0.015 

0 265" l ” 0 ' uu4 
u.zoo _ 0 026 

0.39 ±0.04 

0.35 ±0.04 

bbody 

L BB (.Lzs/DIq) / Norm 

8.8 ± 0.6 

3.6 ±0.8 

5.3 ± 1.3 

4.0 ±0.6 

bbody 

R bb (km) 

7.6 ± 0.5 

17.6 ±4.4 

9.8 ±2.3 

10.6 ±2.5 

nthComp 

kT seed (keV) 

0.90 ± 0.04 

0.569 ±0.014 

0.69 ±0.06 

0.64 ±0.04 

nthComp 

r 

2.05 ± 0.03 

2.081 ±0.018 

2.05 ±0.04 

2.08 ± 0.03 

nthComp 

kT e (keV) 

ioi + 100 

1U -74 

63.2+^ 

80.7+ 5 9 9 9 

89+ 28 

oy -20 

nthComp 

Flux (10 -10 ergs cm -2 s" 1 ) 

4.73 

4.67 

4.66 

4.61 

pexriv 

£ (erg cm s — - 1 ) 

- 

< 1 

< 1 

< 1 

pexriv 

Incl (deg) 

- 

40 (fixed) 

40 (fixed) 

40 (fixed) 

pexriv 

Flux (10 -10 ergs cm " 2 s -1 ) 


0.9 

0.56 

0.73 

gauss 

Dine (keV) 

- 

- 

3.00 ± 0.04 

3.02 ± 0.04 

gauss 

mine (keV) 

- 

- 

0.28 ±0.03 

- 

gauss 

Fine (xio -4 ph cm -2 s _1 ) 

- 

- 

3.3 ±0.9 

2.4 ±0.5 

gauss 

EqW (eV) 

- 

- 

24.0 ±7.6 

18.3 ±3.5 

gauss 

Epe (keV) 

- 

- 

6.21 ±0.08 

6.27 ±0.06 

gauss 

a Fe (keV) 

- 

- 

0.46 ± 0.07 

- 

gauss 

Fine (xlO -4 ph cm -2 s _1 ) 

- 

- 

l-58to.?9 

i Q + 0 - 7 

1>y -0.4 

gauss 

EqW (eV) 

- 

- 

33 ± 17 

52 ±21 

rdblur 

Betor 

- 

- 

- 

-2.2 ±0.5 

rdblur 

Rin ( GM/c 2 ) 

- 

- 

- 

< 29 

rdblur 

Rout (GM/c 2 ) 

- 

- 

- 

360t?®° 

rdblur 

Incl (deg) 

- 

- 

- 

54t ^ 7 

total 

Flux (10 ~ 10 ergs cm -2 s _1 ) 

3.3 ± 0.7 

3.34 ±0.02 

3.34 ± 0.03 

3.34 ±0.03 

total 

X 2 (dof) 

2425 (1511) 

1990 (1509) 

1818 (1503) 

1783 (1501) 


of the fit, since we get y 2 /dof = 1803/1503 fixing the emissivity index at —2.4, the outer 
radius at 400 R 9 , the inclination angle at 40 deg. Also in this case the centroid energy of 
the line remained at 6.1 — 6.3 keV, but we get an upper limit to the inner radius of the disk 
of 25 R g . On the other hand, we get a more significant improvement of the fit if we add 
a (mild) relativistic smearing to the whole reflection component (i.e. the emission lines at 
~ 3 and ~ 6 keV and the pexriv component) convolving all these three components with 
rdblur, the kernel of diskline. In this case we get y 2 /do/ = 1783/1501 (Ay 2 = 35 for the 
addition of three parameters, corresponding to an F-test probability of chance improvement 
of ~ 4 x 10 -7 ). The results of these phenomenological models are reported in Table 1. 

In order to fit the high-quality Suzaku spectrum with more consistent physical models, 
we substit ute th e pex riv + Fe line components with the self-consistent reflection model 
reflionx flRoss & Fabianl 20051). modified by a relativistic blurring component (again mod¬ 
eled with rdblur) to consider the relativistic and/or Doppler effects produced by the motion 
in the inner disk close to the compact object. In this model, emission lines from the most 
abundant elements or ions are also self-consistently calculated. In the reflionx model the 
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Table 2. The best fit parameters of the spectral fitting of the Suzaku spectrum of 4U 1705-44 in the 0.7 — 200 keV energy 
range with the self-consistent reflection model reflionx. The blackbody luminosity is given in units of L 3 $/D 2 q, where L 35 is 
the bolometric luminosity in units of 10 35 ergs/s and Dio the distance to the source in units of 10 kpc. The blackbody radius 
is calculated in the hypothesis of spherical emission and for a distance of 7.4 kpc. The disk blackbody normalization is given 
by (Rin(km) / Dio) 2 cos z, where i is the inclination angle of the binary system. The disk blackbody inner radius R i n (km) is 
calculated for an inclination angle of 40°. Flux is calculated in the 1 — 10 keV band. Uncertainties are given at 90% confidence 
level. 


Parameter 


BBODY [Fe = 1] 

DISKBB [Fe = 1] 

BBODY [Fe = 2] 

DISKBB [Fe = 2] 

Nh (xlO 22 cm - ' 

2 ) 

2.11 ±0.04 

2.27 ±0.04 

2.09 ± 0.04 

2.25 ±0.04 

kT BB (keV) 


0.38 ±0.03 

0.53 ±0.07 

0.38 ± 0.03 

0.52 ±0.08 

L BB {L 35 /D\ o) 

/ Norm 

5.5 ±0.7 

58 ±22 

4.9 ± 0.5 

6 i + 42 

D1 -18 

Rflfl (km) / R irl 

(km) 

10.6 ± 1.8 

6.4 ± 1.2 

10.0 ±1.7 

6 6" 1 ” 2 ' 3 
uo - 1.0 

kT seed (keV) 


0.68 ±0.03 

0.70 ±0.06 

0.66 ± 0.03 

0.67 ±0.06 

r 


2.01 ± 0.02 

2.01 ± 0.02 

1.95 ± 0.01 

1.95 ±0.01 

kT e (keV) 


47+ 19 

47 + 17 
**-10 

27 ±5 

27 ±5 

£ (erg cm s — 1 ) 


< 13 

< 19 

< 24 

22 ±4 

Betor 


-2.5 ±0.5 

-2.5 ± 0.5 

-2.4 ±0.4 

-2.5 ±0.5 

R in ( GM/c 2 ) 


17+ 4 

1 * -6 

16 tj 

16 tt 

16b) 

Rout (GM/c 2 ) 


o'7’/-j-|-8000 

°‘ u -180 

370+ 370 

°‘ u -170 

oonT8000 

oou -200 

35f )”*” 350 

oou_i 70 

Incl (deg) 


43 ±5 

42 ±5 

41 ±4 

42 ±4 

Eline (keV) 


3.03 ±0.03 

3.03 ± 0.03 

3.04 ± 0.03 

3.04 ± 0.03 

hine (XlO -4 ph 

cm -2 s -1 ) 

2.4 ±0.4 

2.3 ±0.4 

2.4 ±0.4 

2.3 ±0.4 

EqW (eV) 


18.2 ±3.6 

17.0 ±3.4 

18.6 ±3.5 

17.5 ±3.7 

Flux (10 -10 ergs 
total x 2 (dof) 

cm -2 s -1 ) 

3.34 ± 0.04 

1831 (1503) 

3.34 ±0.15 

1841 (1503) 

3.34 ± 0.04 

1846 (1503) 

3.34 ±0.14 

1852 (1503) 


emergent (reflected) spectrum is calculated for an optically-thick atmosphere (such as the 
surface of an accretion disk) of constant density illuminated by radiation with a power-law 
spectrum, whose photon index is fixed to that of the nthComp component, and a high-energy 
exponential cutoff with e-folding energy fixed at 300 keV. In order to take into account the 
high-energy cutoff in the illuminating spectrum, we have multiplied the reflionx spectrum 
by a high energy cutoff with the e-folding energy fixed to the value of the e-folding energy 
of the primary (nthComp) component. Since the reflionx model does not take into account 
transitions from Ar and Ca, we fitted the emission line at ~ 3 keV with a diskline, fixing all 
the smearing parameters to those used for the reflionx component. The continuum emission 
is fitted with the same model as before, and we used for the soft component alternatively the 
bbody or the diskbb model. The results of these fits are shown in Table 2. We also checked 


the possib ili 


Egron et a l. 


v of iron overabundance with respect to cosmic abundances (as claimed by 


201 31. for the soft state), fixing the iron abundance alternatively to the cosmic 


value and twice the cosmic value. We find a slightly better fit when we fix the iron abun¬ 
dance to the cosmic value and we use a blackbody to fit the soft thermal component (see 
Table 2). The total 0.5 — 200 keV luminosity of the source during the Suzaku observation 
was 6.15 x 10 36 ergs/s assuming a distance to the source of 7.4 kpc. 

Finally, we fitted the reflection component in the Suzaku spectrum of 4U 1705-44 with the 
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data and folded model 


Unfolded Spectrum 



Figure 3. Left: Suzaku data in the energy range 0.7 - 200 keV (top) and residuals in units of cr with respect to the 
best-fit model (bottom) of 4U 1705-44 (see Table 2, first column). Right: Suzaku unfolded spectrum in the energy range 
0.7 - 200 keV with respect to the best-fit model shown in the first column of Table 2. The model components are also 
shown. From the left to the right we see the blackbody component, the emission line at ~ 3 keV (smeared with the same 
smearing parameters used for the reflection component), the smeared reflection component modeled by reflionx. The main 
Comptonization component and the total model are plotted on top of the data. 


convolution model rfxconv (jKolehma i nen et all 20111), which has the advantage to take into 


account the exact shape of the illuminating continuum. For sake of completeness we also tried 
the self-consistent reflection model relxill by Garcia et_al. (2 0141 ). whose novelty is that 


for each point on the disk the proper reflection spectrum is chosen for each relativistically 
calculated emission angle. In the relxill model we fixed the photon index and e-folding energy 
of the reflected spectrum to those of the nthComp component. Both these models include 
emission lines from the most abundant elements or ions, as in the case of reflionx. Again 
we obtain a good fit of the 4U 1705-44 spectrum, with best-fit values of the reflection and 
smearing parameters well in agreement with those obtained with other reflection models (see 
Tab. 3, cf. with Tab. 1 and 2). 


4 DISCUSSION 


Similar to BH X-ray Binaries, NS LMXBs show clear differences in the spectral parameters 
during hard and soft states. Studying these differences is important in order to address 
different geometry or physical properties of the inner central emitting region in these two 
spectral states and in order to understand what causes the spectral transition. In fact, there 
is no general consensus on the hot corona-accretion disk geometry in these two spectral 
states. Theoretically the hard Comptonized component m ay co me from either a patchy 


corona, possibly powered by m agnet ic flares (e.g. 
located jet (e.g. 


Beloborodov 


Markoff et ah 


19991 ) or the base of a centrally 


2005 ). In both the cases, the thin accretion disk may extend 
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Table 3. The best fit parameters of the spectral fitting of the Suzaku spectrum of 4U 1705-44 in the 0.7 — 200 keV energy 
range with the self-consistent reflection models rfxconv and relxill. The blackbody luminosity is given in units of L 35 /D 1 Q, 
where L 35 is the bolometric luminosity in units of 10 35 ergs/s and D\q the distance to the source in units of 10 kpc. The 
blackbody radius is calculated in the hypothesis of spherical emission and for a distance of 7.4 kpc. Flux is calculated in the 
0.7 — 200 keV band. Uncertainties are given at 90% confidence level. 


Parameter 

RFXCONV [Fe = 1] 

RELXILL [Fe = 1] 

Nfi (xlO 22 cm -2 ) 

2.05 ±0.04 

2.06 ± 0.05 

kTgg (keV) 

0.37 ± 0.03 

0.36 ± 0.04 

L gg ( L 35 /Dj 0 ) / Norm 

4.6 ±0.6 

4.1 ± 0.7 

R bb (km) / R in (km) 

10.2 ± 1.8 

10.3 ± 2.5 

kT seed (keV) 

0.67 ±0.03 

0.63 ± 0.04 

r 

2.01 ± 0.01 

1.979 ± 0.009 

kT e (keV) 

34t® 

43 ±9 

£ (erg cm s — 1 ) 

60l 2 ° 

12 +1 

Refl Amplitude 

0.34 ± 0.04 

0.34 (fixed) 

Refl Norm 

- 

0.161 ± 0.017 

Betor 

-2.5 ±0.5 

o O+0.4 
-0.2 

Rin (GM/c 2 ) 

ntr 

14 ±2 

Rout ( GM/c 2 ) 

> 202 

260 (fixed) 

Incl (deg) 

43 ±5 

31-6 tH 

Eline (keV) 

3.04 ± 0.04 

3.10 ±0.03 

Iline (xio -4 ph cm -2 s -1 ) 

2.5 ±0.4 

2.2 ± 0.4 

EqW (eV) 

19.3 ±0.3 

17±3 

Flux (10 ~ 10 ergs cm -2 s _1 ) 

9.0 ±0.3 

9 33+ 017 

y.oo_0 36 

total x 2 (dof) 

1832 (1503) 

1818 (1504) 


close to the last marginally stable orbit or the NS surface. Alternatively the thin disk may 
be truncated at large distances from the compact object, with the central region replaced 
by an advection-dominated accretion flow (ADAF) where, at high accretion rates, material 
may condense to form an inner optically thick disk (see e.g. Esin et ah 1997). This hot, inner 


flow can also act as the launching site of the jet. In this respect, much information may come 
from the study of the so-called reflection component in different spectral states of a source 


(see 


Done et al. 


20071. as a review) 


In this paper we have reported the results of the spectral analysis of a long Suzaku 
observation of the LMXB of the atoll class 4U 1705-44. This was the result of a ToO 
program intended to observe the source during a hard state. Suzaku observed the source 
for a total of 250 ks yielding a net exposure time of 100 ks. During the observation 9 type- 
I bursts were observed. We present here the spectral analysis of the persistent emission, 
while we will discuss the characteristics of the observed type-I bursts elsewhere. We have 
fitted the persistent emission spectrum in the broad band range between 0.7 and 200 keV 
using the continuum model which gave the best fit of previous high-quality sp ectra o f this 


so urce obtained with RXTE, BeppoSAX, XMM-Newton , and Chandra (see e.g. 


2013 


Egron et al. 


, and references therein). The continuum model consists of a soft component modeled 
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by bbody and a Comptonization component modeled by nthComp, both multiplied by the 
phabs component which takes into account photoelectric absorption by neutral matter in the 
interstellar medium. A smeared reflection component was necessary to obtain an acceptable 
fit of the broad-band Suzaku spectrum. This component was necessary to fit high-energy 
residuals above about 15 keV and a broad iron edge at about 8 — 9 keV - the addition of the 
pexriv component to fit these two features gave a significant improvement of the quality of 
the fit with a probability of chance improvement, calculated with an F-test, which resulted 
to be negligibly small. 

We detect other reflection features, such as a broad emission line at about 6.4 keV from 
neutral iron and a broad emission line at about 3 keV that we tentatively identify with 
the K a transition from Ar XVI-XVII. These two features could be well fitted by broad 
Gaussians. Note that this Ar line is quite strong with respect to the observed iron line. 
If we consider the product of the element abundances (~ 4.7 x 10 5 and ~ 3.6 x 10 6 
for iron and argon, respectively) by the fluorescence yields (which can be calculated using 
the empirical formula yield = Z 4 /(30 4 + Z 4 ), where Z is the atomic number), we obtain 
~ 1.7 x 10 -5 and ~4 x 10 -7 for iron and argon, respectively. This means that the Ar line 
strength should be ~ 2.4% of Fe line strength. This in the hypothesis that all the atoms of 
these elements are in the ionization state producing the line. Note also that the observed 
line may depend on the illuminating continuum at that energy. In the case of 4U 1705-44 
in the soft state, where Ar and Ca lines are clearly detected by XMM-Newton together with 
an Fe line, we find small differences in the line intensities (a factor 2 at most) and higher 
differences in the line equivalent width (up to a factor 8). However, when the spectrum is 


fitted to a self-consisten 


refl ection model, a suitably modified version of the xiliver model 


by Garcia. & Kallman ( 2010 ) which includes Ar and Ca transitions, we find that the Ar 


line is well fitted by the reflection m ode l wit 


i a slight overabundance by a factor 1.8 with 


respect to Solar abundances (see Egron et ah 


2013 ). We conclude therefore that the simple 


calculation above is merely an order of magnitude estimation and that the consistency of 
the Ar line with other reflection features should be checked using self-consistent reflection 
models including Ar and Ca transitions. 

We also tried to fit all the reflection features (the Compton hump, the iron edge and the 
iron line) with a self-consistent reflection model such as ref lionx, which we have modified 
with a high energy cutoff at the electron temperature of the Comptonizing corona to take 
into account the curvature of the Comptonization spectrum with respect to a simple power- 
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Table 4. Comparison of the best-fit continuum and reflection parameters obtained for the soft state (SS) as observed in the 
60-ks XMM-Newton observation and for the hard state (HS ) ob served by Suzaku (this paper). Continuum parameters for the 
SS observed by XMM-Newton are taken by lEgron et al l (' 201 who use a similar model for the continuum, while smearing 
parameters of the reflection component are taken from lDi Salvo et all d2009h where these parameters are o btained with smalle r 
uncertainties. Lx is the X-ray luminosity extrapolated in the 0.1 — 150 keV range for the SS, as reported bv lEgron et alj d2013tl . 
and in th e 0.5 — 2 00 keV range for th e HS (this work). L^dd is the Eddington luminosity for a 1.4 Mg NS, L^dd = 2.5 X 10 38 
ergs s _1 llvan Paradiis fe McClintockll 199411 . 


Parameter 

SS ( XMM-Newton) 

HS ( Suzaku ) 

Nh (xlO 22 cm -2 ) 

2.08 ± 0.02 

2.11 ±0.04 

bbody kT (keV) 

0.56 ±0.01 

0.38 ±0.03 

bbody L x (L 37 /D 2 0 ) 

2.58 ±0.01 

0.055 ± 0.007 

R bb (km) 

33.3 ±1.2 

10.6 ± 1.8 

nthComp kT see d (keV) 

1.30 ±0.02 

0.68 ±0.03 

nthComp kT e (keV) 

3.0 ±0.1 

O O 

<N i-H 
+ 1 
b- 

£ (erg cm s — *) 

> 500 

< 13 

Betor 

-2.27 ±0.08 

-2.5 ±0.5 

R in ( GM/c 2 ) 

14 ±2 

17+ 4 

1 ' -6 

Rout {GM/c 2 ) 

3300~*~ 1500 

oouu_ 900 

Oiyr»-|-8000 

-180 

Incl (deg) 

39 ±1 

43 ± 5 

Lx/L E dd 

72% 

2.9% 


law (used as illuminating spectrum in the reflionx model). The smearing of the reflection 
component has been taken into account multiplying it by the rdblur component. In this case 
we had to add a diskline to the model to fit the Ar line at 3 keV, since transitions from Ar 
are not taken into account in the reflionx model. Anyway, all the smearing parameters of the 
diskline used to fit the Ar line have been fixed to those used for the reflection component. 
In this case, we used a bbody or a diskbb to fit the soft thermal component and tried 
to vary the iron abundance fixing this value to 1 or 2 times the cosmic abundance. The 
best fit corresponds to a bbody component for the soft thermal emission and to an iron 
abundance of 1 (see Table 2). We also tried to fit the reflection component with the self- 
consistent convolution model rfxconv, which consistently takes into account the curvature 
of the illuminating spectrum caused by the high-energy rollover at the electron temperature 
in the corona producing the primary Comptonized spectrum, or with the relxill model, in 
which the cutoff energy in the reflected spectrum is fixed at the electron temperature of the 
primary Comptonized component, fitted with nthComp. Note that both rfxconv and relxill, 
as well as reflionx, all include Compton broadening effects caused by Compton scattering 
in the surface layers of the accretion disc. The results of these fits are reported in Table 3 
and are perfectly consistent with those obtained with all the other reflection models that we 
tried. 

In order to check the stability of the best fit model with respect to the smearing pa¬ 
rameters of the reflection component, we have let all these parameters free to vary. The 
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most uncertain of these parameters, as expected, is the outer radius of the disk, for which 
we find only loose constraints. Interestingly, all of the smearing parameters of the best fit 


model of this observation are in good agreement with the smearing 


obtained with other instruments, e.g. XMM-Newton and BeppoSAX flDi Salvo et al. 


parameters a 


r eady 


JEgron et al 


2009; 


20131 ). during a soft state. In order to facilitate the comparison, we report in 


Table 4 the best fit parameters of the reflection (modeled by reflionx) and the relativistic 
smearing component s ob tained with XMM-Newton during a soft state (from 


Di Salvo et al. 


20091: 


Egron et al.l 2013) and obtained with Suzaku during a hard state (this paper). Al¬ 


though the uncertainties in the spectral parameters in the hard state are larger than in the 
soft state (because of the lower source flux in the hard state), we find a very good agreement 
in all the parameters. In particular, the inclination angles of the system we obtain in the 
two cases are compatible well within the 90% c.l. uncertainty. The main difference in the 
reflection component between the hard and the soft state is in the ionization parameter £, 
which is much larger in the soft state than in the hard state, as expected because of the 
higher incident flux in the soft state. Also the continuum parameters are different; in partic¬ 
ular, the temperatures of the soft components of the continuum (i.e. of the soft blackbody 
component and of the seed photons for the Comptonization) result higher in the soft state, 
while the electron temperature of the Comptonizing cloud results higher in the hard state, 
in agreement with what is expected. 

We note that there is no clear indication of a receding inner accretion disk radius in the 
hard state, corresponding to a luminosity of ~ 3% of the Eddington luminosity, with respect 
to the soft state, which was observed at about 70% of L Edd . On the contrary, the inner 
disk radius as inferred from the reflection component is consistent to be the same in the two 
spectral states, at about 34 km from the NS center. A similar indication comes from the inner 
radius of the disk as inferred from the blackbody component, that we interpret as the direct 
emission from the accretion disk. Both in the soft and in the hard state the blackbody radius 
is a few tens of km, in agreement with the estimate we get from the reflection component. We 
caution the reader, however, that neither the color factor or the geometry of the system has 
been taken into account in this calculation. What is reported is just a zero-order estimation 
of the radius of the region (assumed to be spherical) emitting the blackbody component. In 
particular, the spectral hardening factor may depend on luminosity (see e.g. 


Merloni et al. 


2000) explaining why the inner disk radius may appear l arger at 


ligher luminosities. 


This result is in agreement with what is found by lEg ron et al. (12013 1 who studied XMM- 
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Newton, BeppoSAX, and RXTE spectra of 4U 1705-44 in the hard state and in the soft 
state. In particular, in the hard state, the inner disk radius as measured by the smearing 
of the reflection component resulted at 19 — 59 K g , which is compatible with the inner 
disk radius derived in the soft state (13 ± 3 R s ), while more uncertain results came from 


the evaluation of the blackbody radius in the hard state. 


D’Ai ct al 


the same XMM-Newton observation during a hard state used bv lEgron et ah! (120131 ). These 


(2010h also a nalyzed 


authors discussed the possibility of a truncated disk in the hard state based on the best fit 
value of the inner disk radius as found from the Fe line width, which was about 30 R g , that 
means about 60 km. However, in that case, the lack of broad-band coverage and the limited 
statistics, gave a large uncertainty on the inner disk radius, whose 90% c.l. range was from 
6 to 90 R r Considering the large uncertainty on this measurement we cannot state that 


the result was in contrast with more recent results ( Egron et al 


also that the best-fit blackbody radius reported by 


D’Ai ct al. 


2013, and this work). Note 


(120101) has a value around 


14 ±5 km , indicating that the disk may be truncated quite close to the NS surface. Similarly, 
( 120101 ) could not determine with high precision the inner radius of the disk using a 


Lin et al. 


diskline model for the Fe line in Suzaku and BeppoSAX spectra of 4U 1705-44 taken during 
a hard state. Therefore, the inner radius of the disk was fixed to 6 R s , and the fit results, 
such as Fe line flux and equivalent width, were not sensitive to this parameter. 

Similar results for the inner disk radius were obtained also in the case of 4U 1728-34, 


the prototype of the atoll sources. The XMM-Newton spectrum reported by [Egron et al 


(2 0111) taken during a low-luminosity state of the source (probably a hard state) showed a 
relatively broad iron line (Gaussian a ~ 0.6 keV), which was fitted to a series of models 
(diskline, relline, and reflionx) yielding in all the cases an inner disk radius b etwee n 12 and 


22 R„. In this case a blackbody component was not significantly detected. 


Cackett et al. 


( 2010 ) present a comprehensive, systematic analysis of Suzaku and XMM-Newton spectra 
of 10 NS LMXBs, in order to study their Fe Ka emission lines. In most cases they find a 
narrow range of inner disk radii (6 — 15 R g ), implying that the accretion disk extends close 
to the NS surface over a range of luminosities. 

In this respect, it may be useful to compare these results to those obtained for BH X- 
ray Binaries, since much work has been done to determine the inner radius of the disk in 
these systems bot 
component (e.g. 


:h from the iron 

ine and the reflection component and from the blackbody 

Done et al. 

2007; 

Reis et al. 

2010 

, and references therein). Also for these 

growing evidence that the disk may not be truncated far from the last 
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stable orbit. The broad-band (0.1 — 200 keV) BeppoSAX spectrum of one of the best studied 
galactic BH candidates, Cygnus X-l, taken during a hard state showed evidence of a complex 


reflection component. In this spectral deconvolution the the inner rad ius o f the c 


inferred from the smeared reflection, is found between 6 and 20 R,„ (IDi Salvo et al. 


isk, as 


This result is in agreement with the results of 


Young et ah 


200111 . 


(2001) who fitted ASCA, Ginga 


and EXOSAT data of Cygnus X-l in both soft and hard spectral states to a model of 
an ionized accretion disk, whose spectrum is blurred by relativistic effects. They found 
that relativistic blurring provided a much better fit to the low/hard state and that data 
of both states were consistent with an ionized thin accretion disk with a reflected fraction 
of unity ext ending to the innermost stable circular orbit around the BH (see, however, 
for a different interpretation). Up to date, one of the strongest evidence 


Barrio et ah 


2003 


of a trunc ated disk, base d on Fe line profile measurements, in a BH hard state is that of 
GX 339-4 (ITomsick et ah 20091). In that case, from Suzaku and RXTE spectra, it was found 
that R m was a factor > 27 higher than in the bright state when the luminosity was about 
at 0.14% of the co rresp onding Eddington limit. 

More recently Reis et ah (2010) have analyzed a sample of stellar mass BHs, including 
Cygnus X-l, in the low-hard state, down to luminosities of ~ 10~ 3 L Edd , finding no clear 


evidence of a truncation of the inner disk at radii larger than 10 R 3 . Furthermore, the 
thermal-disk continuum yields colour temperatures consistent with the relation L oc T 4 , 
implying that the emitting surface is consistent with being constant with luminosity. A 
similar relation, L oc T 3,2 , seems also to hold in the case of 4U 1705-44 at least for the 


soft state (Lin et ah 


2010). The authors suggest that the deviation may be caused by a 


luminosity-dependent spectral hardening factor. On the other hand, the relatively low re¬ 
flection amplitude (0/27 t ~ 0.34) we find in the hard state of 4U 1705-44 fitting its spectrum 
to self-consistent reflection models, such as rfxconv, is compatible with a spherical geometry 
with the hot (spherical) corona inside an outer accretion disk. In this case, the small inner 
radius of the disk we find (R in ~ 17 R s , corresponding to approximately 30 — 35 km for a 1.4 
M 0 NS) would indicate a very compact hot corona filling the central part of the accretion 
disk, which in the case of a NS may be identified with a boundary layer between the inner 
accretion disk and the NS surface. This might represent an important difference between 
BH and NS systems, since for a BH, in the absence of a boundary layer, the inner disk 
should extend down to the last stable orbit. Note, however, that alternative explanations 
for this weak reflection cannot be ruled out. This can be caused by a highly ionized inner 
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disk 

disk 

e.g. 

Ross et al. 

1999 

) or mildly relativistic outflow of the hot corona away from the 

Beloborodov 

1999 

). In the latter case, it was shown that reflection fractions as low as 


~ 0.3 can be obtained in the low-hard state without invoking a truncated disk. 


In summary, we have analyzed a deep (100 ks exposure time) ToO Suzaku observation 
of 4U 1705-44 during a low-luminosity hard state (corresponding to a luminosity of ~ 3% 
L Edd)- The broad-band spectrum shows a prominent Compton hump at hard energies, a Fe 
absorption edge and two relatively weak emission lines at ~ 3 and ~ 6.4 keV, identified 
with fluorescent emission from mildly ionized Ar and neutral Fe, respectively. We used all 
the available to date self-consistent reflection models to fit these reflection features finding 
in all cases best-fit parameters that are compatible with each other and consistent with 
those reported in literature for the soft state. In particular the inclination angles found from 
the reflection component for the hard and the soft state are perfectly compatible with each 
other. It is worth noting that we obtain similar smearing parameters in the soft and in the 
hard state even if in the soft state the reflection component is dominated by the Fe line, 
while in the hard state other features are dominant. In fact the Fe line is the least important 
feature in the statistical sense (since the addition of other reflection features, such as the 
Compton hump plus the Fe edge, gave the most important improvement of the fit). We also 
find very similar inner disk radii in the hard and soft state, indicating that the inner disk 
rim does not change significantly at different spectral states down to a luminosity of ~ 3% 
the Eddington limit. 
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